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Abstract - Recent collections of specimens from the Pilbara attributed to 
Leggadina lakedownensis stimulated an investigation of genetic and 
morphometric variation in L. lakedownensis and L.forresti to confirm the 
taxonomic status of these species, determine the extent of intraspecific 
variation, and develop useful markers for classification. Extensive analyses of 
both morphology and allozymes confirm the clear taxonomic separation of 
these two species. Leggadina lakedownensis shows considerably more 
intraspecific morphological variation than LJorresti, reflecting the former's 
much greater geographic range and ecological diversity. The mean body size 
of Thevenard Island L. lakedownensis is markedly larger than the other 
populations examined. Despite this morphological diversity, populations of 
L lakedownensis, including those from mainland Pilbara and Thevenard 
Island, are genetically similar across their range. An individual from 
Queensland is genetically the most divergent but we have been unable to 
assess the morphological relationship of the topotypical population. 


INTRODUCTION 

The genus Leggadina has a wide distribution 
across mainland Australia, with two allopatric 
species currently recognised, L.forresti (Thomas, 
1906) from central Australia and L. lakedozvnensis 
(Watts, 1976), broadly distributed across northern 
Australia (Figure 1). Baverstock et al. (1976) 
presented biochemical and karyotypic evidence to 
firmly establish the specific distinction of 
L. lakedozvnensis, then known only from a restricted 
area of northeast Queensland, from the more widely 
distributed Leggadina forresti. 

Since these initial studies, many different 
populations of Leggadina have been sampled, 
especially in the Kimberley and Pilbara regions of 
Western Australia. Although some of these 
populations have been referred to L lakedozvnensis, 
this attribution has not been critically examined. In 
addition, animals from Thevenard Island, off the 
Pilbara coast, are much larger than mainland 
populations of Leggadina, and this has led to some 
concern over the taxonomic status and conservation 
needs of this morphologically-distinctive form. 
Moro et al. (1998) examined variation in 


mitochondrial DNA (mtDNA) from Western 
Australia and Northern Territory populations 
referred to Leggadina lakedozvnensis and found 
significant lineage divergence between the Pilbara 
(including Thevenard Island) and Kimberley 
populations of this species. On this basis they 
suggested that the Pilbara and Thevenard Island 
specimens together might represent a distinct 
subspecies. Moro et al. (1998) did not include any 
specimens from northeast Queensland in their 
analyses. 

Prior to 1997, only two specimens of Leggadina 
had been collected in the Pilbara region. Since that 
time, thirty specimens, some with tissue for 
molecular analysis, have been collected during 
various surveys of the region. The availability of 
this new material has prompted us to conduct a 
comprehensive morphological and genetic study of 
Leggadina across its geographic range. The purposes 
of this study were to re-assess species boundaries in 
the genus using morphological and allozyme data, 
and to examine broader patterns of morphological 
and genetic variation across the range of L. 
lakedozvnensis. 
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Figure 1 Distribution of Leggadina, with localities for specimens used in this study marked, (A) L lakedownensis, 
(•) L.forresti. 


MATERIALS AND METHODS 

Allozyme characterization 
Genetic analysis was based on tissue samples 
from 16 specimens of Leggadina lakedownensis, 
aggregated into six populations according to 
geographic location, and 12 specimens of L.forresti 
(Appendix 1). Liver samples for allozyme 
electrophoresis were removed from freshly killed 
animals in the field or at the Western Australian 
Museum and immediately frozen in liquid nitrogen 
before being stored at -SO^C. Allozyme 
electrophoresis of liver homogenates was 
conducted on cellulose acetate gels ("Cellogel') 
according to the methods of Richardson et al. (1986). 
The following proteins exhibited zymograms of 
sufficient activity to assign genotypes: aconitase 
hydratase (ACON, EC 4.2.1.3), aminoacylase 
(ACYC, EC 3.5.1.14), adenosine deaminase (ADA, 
EC 3.5.4.4), alcohol dehydrogenase (ADH, EC 
1.1.1.1), adenylate kinase (AK, EC 2.7.4.3), albumen 
(ALB), carbonate dehydratase (CA, EC 4.2.1.1), 
diaphorase (DIA, EC 1.6.99.), enolase (ENOL, EC 
4.2.1.11), fructose-bisphosphatase (FDPASE, EC 


3.1.3.11), fumarate hydratase (FUM, EC 4.2.1.2), 
glyceraldehyde-3-phosphate dehydrogenase 
(GAPD, EC 1.2.1.12), lactoylglutathione lyase (GLO, 
EC 4.4.1.5), aspartate aminotransferase (GOT, EC 
2.6.1.1), glucose-6-phosphate dehydrogenase 
(G6PD, EC 1.1.1.49), glycerol-3-phosphate 
dehyrogenase (GPD, EC 1.1.1.8), glucose-6- 
phosphate isomerase (GPI, EC 5.3.1.9), glutathione 
peroxidase (GPX, EC 1.11.1.9), glutathione 
reductase (GSR, EC 1.6.4.2), isocitrate 
dehydrogenase (IDH, EC 1.1.1.42), L-lactate 
dehydrogenase (LDH, EC 1.1.1.27), malate 
dehydrogenase (MDH, EC 1.1.1.37), '"malic" 
enzyme (ME, EC 1.1.1.40), mannose-6-phosphate 
isomerase (MPI, EC 5.3.1.8), purine-nucleoside 
phosphorylase (NP, EC 2.4.2.1), dipeptidase (PEP- 
A, EC 3.4.13.), tripeptide aminopeptidase (PEP-B, 
EC 3.4.11.), proline dipeptidase (PEP-D, EC 3.4.13.), 
phosphoglycerate mutase (PGAM, EC 5.4.2.1), 
phosphogluconate dehydrogenase (6PGD, EC 
1.1.1.44), phosphoglycerate kinase (PGK, EC 
2.7.2.3), phosphoglucomutase (PGM, EC 5.4.2.2), 
pyruvate kinase (PK, EC 2.7.1.40), superoxide 
dismutase (SOD, EC 1.15.1.1), and triose-phosphate 
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isomerase (TPI, EC 5.3.1.1). The nomenclature for 
designating allozymes and multiple loci follows 
Adams et al. (1987). 

A phylogenetic analysis using all specimens from 
both species was undertaken by constructing a 
neighbor-joining tree (Saitou and Nei, 1987) of Nei 
(1978) unbiased genetic distance. Because the 
populations of L. lakedownensis could potentially be 
interbreeding we also investigated their genetic 
relationships by a principal coordinates analysis 
(Gower, 1966) using the Cavalli-Sforza and 
Edwards (1967) arc distance matrix which is 
Euclidean. These analyses were carried out using 
Genstat (2000), Biosys-1 (Swofford and Selander, 
1989) and Mega version 2.1 (Kumar et aL, 2001). 

Morphological characterization 

Morphological analysis was based on 132 
specimens, of which 108 were attributable to 
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L. lakedoxvnensis and aggregated into eight 
geographic populations, and the remainder (24) 
spread across most of the geographic range of 
L forresti (Appendix 2). Skulls of L lakedownensis 
from Queensland were not available for 
measurement until after the morphological analysis 
was completed. From the analysis, the significance 
of the length of the upper molar row was realised 
and then these measurements of the Queensland 
specimens were taken, though this group was not 
used in any multivariate statistical analyses. Skulls 
from Fowlers Gap and Innamincka (SAM13633, 
SAMI0170) identified as L forresti by Baverstock et 
al, (1976) using karyotypic and biochemical analyses 
were included herein, as were skulls from 
Oodnadatta, Birdsville, MacDonald Ranges and 
Southwest Mann Range (M21623, M21762, M21786, 
M6344, and M6345), also identified as L forresti by 
Watts (1976) using morphology. 





". BL 

Figure 2 Skull measurement points. 
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Twenty measurements of skull, dentary and 
dental characters and six measurements of external 
characters (Figure 2) were recorded to the nearest 
0.01mm using digital callipers. The measurements 
were GL, greatest skull length; BD, braincase depth; 
IPW, interparietal width; BW, braincase width; ZW, 
zygomatic width; IW, interorbital width; NL, nasal 
length; NW, nasal width; MW, mastoid width; 
BUW, width across bullae; BL, basal skull length; 
M^L, upper molar 1 length; M^L, upper molar 3 
length; MW, upper molar 1 width; M^M^L, length 
of upper molar row; APF, anterior palatal foramen 
length; PL, palatal length; CIL, dentary condyle to 
tip of incisor; CAL, tip of coronoid process to 
posterior edge of angular process; M^M^L, lower 
molar row length. The external measurements were 
HV, head to vent; TV, vent to tail tip; TAILW, tail 
diameter measured close to body; PESW, width of 
pes at hallux; PESL, pes length; EARL, ear length. 

M‘L was not measured during this study because 
its significance in differentiating between 
L lakedozvnensis and Lforresti was not appreciated 
at the time of measurement, it was later calculated 
as M'M'L-M'L-M^L. 

Individuals were categorised into three age 
groups, mature, adult and juvenile according to the 
extent of tooth-wear and fusion of cranial sutures; 
mature ~ extensive tooth wear and sutures fused, 
adult- some tooth wear and most sutures fused and 
juveniles- tooth eruption, no tootli wear and sutures 
unfused. Age was confirmed by examination of the 
external reproductive status of the animal, fully 
descended testes or nipple development 
representing adults and minimal testes descent and 
button nipples representing juveniles. Three 
individuals were classified into juvenile and were 
excluded from all analyses. Appendix 3 lists 
summary statistics for 20 skull and 6 external 
variables in L. lakedozvnensis (grouped by 
geographical region) and Lforresti, (all localities 
pooled). 

Initially, we used multiple regressions to 
investigate and remov'e any effects of age and sex 
before further multivariate analyses. These analyses 
also permitted the detection of statistical outliers 
and non-normality. For L. lakedozvnensis, each 
character was regressed on sex, age and locality, 
and the interactions between these factors. The 
interactions were examined first because they test 
for the inter-dependence of the three factors. CAL 
was the only cranial variable that had a significant 
interaction (between age and locality). This was a 
consequence of mean differences between the adults 
and mature age classes varying between localities. 
There was also a significant two-way interaction for 
TAILW (p < 0.01); there was no obvious fat storage 
in tails. These two variables were not used in any 
further multivariate analyses. Four variables 
showed statistically significant sexual dimorphism 
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at the 5% level of significance (see results) and most 
variables revealed statistically significant 
differences between age classes and between 
localities. We therefore used the estimates of the 
regression coefficients to remove the age and sex 
differences while retaining the inter-locality 
variation. Male and female values were adjusted to 
the mean of the two sexes, and mature and adult 
values to the mean of these two age classes. For 
Lforresti, there were not enough individuals of both 
sexes and each age group to test for interactions or 
reliably estimate main effects, so sexual 
dimorphism and age differences in this species 
were corrected by using the regression estimates 
made from L. lakedozvnensis. All analysis were 
repeated with Lforresti variables corrected for sex 
and age effects using regression estimates made 
from L. forresti and results were unchanged. 

Canonical variate analysis (discriminant analysis) 
was used to determine the multivariate 
relationships between groups. Where a subset of 
characters was required, they were selected by 
sequential multivariate analysis of variance, using 
backward elimination and Wilks' lambda as the 
selection criteria. All statistical analyses were 
carried out using Genstat (2000). 


RESULTS 

Interspecies comparison 

Allelic variation was observed at 24 of the 44 
putative genetic loci examined (Table 1). The 
genetic data unequivocally diagnosed 
L. lakedozvnensis from L. forresti, with fixed 
differences at seven loci {AconJ, Enol, Fdp, G6pd, 
Gotl, Gsr, and Mdh2) and major differences in allele 
frequency at another three loci {Adbl, Gda, and Sod). 
The unbiased Nei Distance between the species is 
about 0.28, indicating marked genetic 
differentiation (Table 2 and Figure 3). In contrast, 
regional populations of L lakedozvnensis showed 
minimal differentiation, with no pairwise Nei 
Distance greater than 0.033 (Table 2), and only one 
locus showed any evidence of genetic heterogeneity 
(ie. a deficiency of heterozygotes) in Lforresti (Ada, 
data not shown). 

The initial morphological appraisal examined the 
relationship between L. lakedozvnensis and Lforresti 
by canonical variate analysis of cranial measures to 
confirm the specific status of these taxa. While there 
is no overlap in the distribution of the two species 
on the canonical variate, two L lakcdozvnezzsis are 
positioned closer to the Lforresti centroid than their 
own species and so are incorrectly allocated to the 
wrong species (Figure 4). Because the sample size 
of L forresti was small, only five traits, selected by 
backward elimination, were used in this analysis. 
The canonical variate is a contrast of upper molar 1 
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Table 1 Sample size for genetic analysis (N), putative allele frequencies at 24 variable loci, and mean heterozygosity 
(H) in 6 populations of L. lakedotvnensis and in Lforresti. The following 20 loci were invariant:- Akl, Akl, Alb, 
Fum, Gapd, Glo, Gpd, Gpi, Gpx, Idhl, Ldhl, Ldh2, Mdhl, Me, Np, PepAl, PepB, Pgk, Pgml, and Tpi 

L. lakedownensis L. forresti 


Locus 

Allele 

Pilbara 

Thevenard 

Island 

Mitchell 

Plateau 

Lissadell 

Kakadu 

Queensland 


N 


6 

1 

3 

2 

3 

1 

12 

Aconl 

a 

- 

- 

- 

- 

- 

- 

1.00 


b 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

- 

Aconl 

a 

- 

- 

- 

- 

- 

- 

0.04 


b 

1.00 

1.00 

0.83 

1.00 

1.00 

1.00 

0.88 


c 

- 

- 

0.17 

- 

- 

- 

0.08 

Acyc 

a 

- 

- 

0.17 

- 

- 

- 

0.04 


b 

0.92 

1.00 

0.83 

0.75 

0.67 

1.00 

0.92 


c 

0.08 

- 

- 

0.25 

0.33 

- 

0.04 

Ada 

a 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

0.42 


b 

_ 

_ 

_ 

_ 

- 

- 

0.50 


c 

_ 

- 

- 

- 

- 

- 

0.08 

Adhl 

a 

- 

- 

- 

- 

- 

- 

0.88 


b 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

0.12 

Adh 2 

a 

0.08 

- 

- 

- 

- 

- 

- 


b 

0.92 

1.00 

1.00 

0.75 

0.83 

1.00 

1.00 


c 

- 

- 

- 

0.25 

0.17 

- 

- 

Ca 

a 

0.17 

- 

- 

- 

- 

- 

0.12 


b 

0.83 

1.00 

1.00 

1.00 

1.00 

1.00 

0.88 

Dia 

a 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

0.92 


b 

- 

- 

- 

- 

- 

- 

0.08 

Enol 

a 

- 

- 

- 

- 

- 

- 

1.00 


b 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

- 

Fdp 

a 

- 

- 

- 

- 

- 

- 

1.00 


b 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

- 

GSpd 

a 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

- 


b 

_ 

- 

- 

- 

- 

- 

1.00 

Gda 

a 

- 

- 

- 

- 

- 

- 

0.88 


b 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

0.12 

Got! 

a 

- 

- 

- 

0.50 

- 

- 

1.00 


b 

1.00 

1.00 

1.00 

0.50 

1.00 

1.00 

- 

Got2 

a 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

0.96 


b 

- 

- 

- 

- 

- 

- 

0.04 

Gsr 

a 

0.08 

- 

- 

- 

- 

- 

- 


b 

_ 

- 

- 

- 

- 

- 

1.00 


c 

0.92 

1.00 

1.00 

1.00 

1.00 

1.00 

- 

ldh2 

a 

0.92 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 


b 

0.08 

- 

- 

- 

- 

- 

- 

Mdh2 

a 

- 

- 

- 

- 

- 

- 

1.00 


b 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

- 

Mpi 

a 

_ 

_ 

- 

- 

- 

- 

0.08 


b 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

0.92 

PepA2 

a 

- 

- 

0.17 

0.25 

- 

- 

- 


b 

1.00 

1.00 

0.83 

0.75 

1.00 

1.00 

1.00 

PepD 

a 

- 

- 

0.33 

- 

- 

- 

- 


b 

1.00 

1.00 

0.67 

1.00 

1.00 

1.00 

1.00 

Pgam 

a 

- 

- 

- 

- 

- 

- 

0.04 


b 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

0.96 

6Pgd 

a 

_ 

- 

- 

0.25 

0.67 

- 

0.75 


b 

1.00 

1.00 

1.00 

0.75 

0.33 

0.50 

0.25 


c 

_ 

- 

- 

- 

- 

0.50 

- 

Pk 

a 

1.00 

1.00 

1.00 

1.00 

1.00 

- 

1.00 


b 

_ 

- 

- 

- 

- 

1.00 

- 

Sod 

a 

- 

- 

- 

0.25 

- 

- 

1.00 


b 

1.00 

1.00 

1.00 

0.75 

1.00 

1.00 

- 

H 


0.022 

0.000 

0.035 

0.072 

0.032 

0.023 

0.058 
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Figure 3 Neighbor-joining tree based on unbiased Nei genetic distance between six populations of L. lakedozvnensis 
and L forresti. 


Table 2 Genetic distances amongst populations of L lakedozvnensis plus L forresti. Nei unbiased distance at bottom left 
and Cavalli-Sforza and Edwards arc distance above and right. 



Pilbara 

Thevenard 

Island 

Mitchell 

Plateau 

Lissadell 

Kakadu 

Queensland 

L. forresti 

Pilbara 

_ 

0.069 

0.115 

0.142 

0.122 

0.182 

0.481 

Thevenard Island 

0.000 

- 

0.092 

0.135 

0.116 

0.169 

0.481 

Mitchell Plateau 

0.002 

0.001 

- 

0.150 

0.147 

0.192 

0.486 

Lissadell 

0.003 

0.003 

0.002 

- 

0.112 

0.215 

0.447 

Kakadu 

0.010 

0.011 

0.011 

0.001 

- 

0.200 

0.473 

Queensland 

0.030 

0.029 

0.031 

0.031 

0.033 

- 

0.507 

L. forresti 

0.278 

0.278 

0.285 

0.236 

0.266 

0.308 

- 


length (M^L) with the other four traits (Table 3). An 
equivalent analysis using five external measures 
revealed poor separation of the two species. 
Furthermore, 20 of 78 L lakedownensis and one of 
eight L. forresti specimens were incorrectly 
allocated. 

Geographic variation within L. lakedownensis 
Principal coordinates analysis of the six 
L. lakedozvnensis populations revealed three genetic 


clusters, namely (a) Queensland, (b) Lissadell and 
Kakadu, and (c) the Pilbara, Thevenard Island and 
Mitchell Plateau populations (Figure 5). This 
outcome can be seen in the genetic distance Vcilues, 
with minimal genetic differentiation between the 
three Western Australian populations and between 
the two Northern Territory populations (all Nei 
Distances < 0.003, Table 2). 

Multiple regressions of skull and external 
measurements rev^ealed that only four variables 


L. lakedownensis 

L. forresti i 

II 1 HIM 

1 1 II 1 1 

1 1 1 

in II Riiiiiiiiiiiiiiii III iim 1 1 

nil II II 1 1 

1 

-€ 

1 

-4 

1 

-2 

Canonical Variate 

0 

1 

2 


Figure 4 Plot of individuals on the canonical variate discriminating L. lakedownensis and L forresti. 
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Figure 5 Scatter plot of first two components from a 
principal coordinates analysis of the arc 
genetic distance matrix for six L lakedownensis 
populations. 


Tables Correlations between original cranial 
measures and canonical variate for 
discriminant analysis of L. lakedownensis and 
L.forresti. 


Trait 

Correlation 

BD 

-B.43 

BL 

-0.35 

M'L 

0.36 

M'M^L 

-0.12 

CIL 

-0.12 


showed statistically significant sexual dimorphism 
(0.01 < p < 0.05 for all four variables). The male 
MW average was 0.16mm less them females; male 
BUW average was 0.15mm less than females; and 
males had a broader and shorter pes compared to 
females (PESW 0.19 larger and PESL 0.48 smaller). 
Most morphometric variables revealed statistically 
significant differences between age classes and 
between localities. All variables except APE showed 
statistically significant differences between 
localities, most at the 0.001 level. Individuals of 
L. lakedownensis from Thevenard Island and Pilbara 
mainland localities were larger than individuals 
from every other locality (Appendix 3). 

Canonical variate amalysis using cramial measures 
for the seven localities revealed the Pilbara and 
Thevenard island populations clustered separately 
from the other five locations in northern Australia 
(Figure 6). This analysis was based on five traits 
selected by backward elimination, but essentially 
the same results were obtained using all 19 cranial 
measures or forward selection of the five best 
(which produced four of the same traits as 


Table 4 Correlations between cranial measures and the 


first and second canonical variates (CV 1 and 
CV 2) for seven L lakedoivnensis localities. 


Trait 

CVl 

CV2 

BUW 

0.71 

-0.33 

BL 

0.53 

0.24 

M’L 

0-30 

0.26 

M’W 

0.74 

0.57 

CIL 

0.38 

0.36 


Table 5 Correlations between external measures and 
the first and second canonical variates (CV 1 
and CV 2) for seven L. lakedownensis localities. 


Trait 

CVl 

CV2 

HV 

0.21 

-B.28 

TV 

0.79 

-0.54 

PESW 

0.38 

0.43 

PESL 

0.66 

0.46 

EARL 

0.13 

0.31 


backward selection). The first two canonical 
variates accounted for 73% and 17% of the between 
group variance (Table 4). The uniformly positive 
correlations between the fiv^e traits and the first 
canonical function show that this is largely a size 
discriminator, accounting for the separation of the 
Pilbara and Thevenard Islamd populations from the 
others on this axis. The second canonical variate is a 
contrast of width across bullae and the other four 
traits, especially width of Mh 
A parallel cinalysis using the five external features 
clearly separated the Thevenard Island population 



1st Canonical Variate 

Figure 6 Plot on the first two canonical variates, using 
skull variables, of seven L. lakedoivnensis 
populations showing group centroids and 
polygons enclosing the range of individals. 
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1st Canonical Variate 

Figure 7 Plot on the first two canonical variates, using 
external variables, of seven L lakedoionensis 
populations showing group centroids and 
polygons enclosing the range of individals. 
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from the others on the first canonical variate (Figure 
7). This axis accounts for 85% of the be tween-group 
variation and is positively correlated with all 
measures (Table 5), this reflects the large size of the 
Thevenard Island animals, particularly in regard to 
the length of the tail-vent and pes. 

Morphological variation in L. lakedownensis and 
L,forresti 

Skull and teeth 

Differences in skull and tooth morphology of 
L lakedoxvnensis and Lforresti are summarised in 
Table 6. Figures 8a-<:, are dorsal, ventral and lateral 
photographs of skulls of L. lakedoxvnensis, the 
smaller skull, M16906 and Lforresti, the larger 
skull, M6352. Both species appear quite variable, 
witli especially marked geographic variation in the 
former species. The molar structure of both species 
is very similar with v^ariation in the relative size of 
individual cusps and teeth rather than in their 
overall shape. From the canonical variate analysis, 
the variables BD, BL, M^L, M'M^L, CIL were the 
most significant in diagnosing the two species; these 



Figure 8a Dorsal view of L. lakedownensis, left, M16906 and Lforrestii, right, M6352. 
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Figure 8b Ventral view of 
L. lakedownensis, 
left, M16906 and 
L.forrestii, right, 
M6352. 



Figure 8c Lateral view of 
L. lakedownensis, below. 
Ml6906 and L.forrestii, 
above, M6352. 


V 
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Table 6 Skull character differences between L lakedownensis and L forresti. 
L. lakcdozvnensis L. forresti 


Upper incisor sometimes proodont 
Rostrum shorter, broader 
Braincase "square" dorsally 
Braincase shallower 

Basicranium shorter (except for Thevenard Island animals) 
Occipital surface of skull relatively flat, 
with weak nuchal crest 
Nasal bone stops short of incisor 
Bullae wider, longer, flatter 
M'L generally longer 
M‘ anterior cingulum long 
shorter relative to 
Mj longer 

Dentary angular process large 
Condyle to tip of incisor generally shorter 


Upper incisor not proodont 
Rostrum longer, narrower 
Braincase "round" dorsally 
Braincase deeper 
Basicranium longer 

Occipital surface of skull with distinct bulge in 

supracoccipital region, nuchal crest well-developed 
Nasal bone usually overlaps incisor 
Bullae shorter, narrower, more inflated 
M'L generally shorter 
M' anterior cingulum short 
longer relative to M' 

Mj shorter 

Dentary angular process finer 
Condyle to tip of incisor generally longer 


are listed as distinguishing characters in Table 6. 
The two species differ iia several cranial and dental 
charcteristics (Tables 6 and 7). One particularly 
diagnostic feature is the relative length of M' to 
(M^ is relatively longer in L. lakedownensis). This 
reflects the relatively greater development of the 
M' anterior cingular cuspule of L. lakedoivnensis, 
which approximates cusp tl in bulk in this species 
but is much smaller than tl in L forresti. This 
difference is reflected in a different overall crown 
outline for the M', with a strongly L-shaped 
anterolingual inflection in L. lakedownensis and a 
more gently concave inflection in L. forresti. 
Variation in cingular morphology is also observed 
within L. lakedownensis. In Kimberley specimens of 
LAakedownensis and in L. forresti the anterior 
cingulum extends posteriorly around the lingual 
face of cusp tl. In contrast, in Pilbara and 
Thevenard Island specimens, the cingulum 
terminates on the anterior face of cusp tl. 

Across the range of L. lakedoivnensis the incisive 
foramina only rarely widen posteriorly (Watts, 


1976), rather thmi narrowmg as is usual in L. forresti 
and there is only rarely a notch on the upper incisor 
(Grahame Medlin, pers comm.). No notches on the 
incisor have been observed in L forresti. 

The skull morphology of L. lakedownensis is 
variable across its geographic range, with the upper 
incisor usually proodontic in Kimberley 
populations and only occasionally in others; the 
Pilbara and Thevenard Island animals have larger 
skull, dental and dentary characters than other 
populations; and the accessory cusp in the 
Thevenard Island animals is not as elongated as in 
other groups. 

Externals measurements 

Mean values of external characters m L forresti are 
larger than those of all L lakedoivnensis populations 
except Thevenard Island, where the animals 
approach the size of L. forresti (Appendix 3). An 
exception is mean ear length which distinguishes 
the two species absolutely. Pes of L. forresti are not 
significantly wider than L lakedownensis, contrary to 


Table 7 Table of means and ratios of M’ and for adult male and female L. forresti and L lakedoivnensis. 




L. forresti 


L. 

lakedownesis 

M^L 

M'L 

%MVM' 

M^L 

M'L 

%MVM' 

Male Total 

1.35 

2.24 

60.3% 

1.19 

2.40 

49.6% 

Female Total 

1.40 

2.20 

63.6% 

1.23 

2.40 

51.3% 

Pilbara male 




1.33 

2.48 

53.6% 

Pilbara female 




1.29 

2.53 

50.1% 

Thevenard male 




1.28 

2.49 

51.4% 

Thevenard female 




1.26 

2.43 

51.9% 

Kimberley male 




1.15 

2.38 

48% 

Kimberley female 




1.21 

2.37 

51% 

Queensland female 




0.82 

2.47 

46% 
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Watts (1976). Mean tail length relative to head and 
body length for L lakedoivnensis ranges from 72% in 
males in the Northern Territory to 91% for males in 
the Pilbara and 99% for females in Thevenard. For 
L.forresti the mean values are 81% for females and 
83% for males. 

Pelage and footpads 

Unlike L.forresti, L. lakedoivnensis is variable in 
pelage colour, with the palest being Queensland 
specimens and the darkest being from the Pilbara. 
L.forresti is more rufous dorsally and it has a pure 
white belly. L. lakedoivnensis is more mouse-like 
(greyish) in colour overall, with only Queensland 
and Great Sandy Desert specimens having a white 
to cream coloured belly, though never as starkly 
white as in L.forresti. The tails of L.forresti are 
generally hairier in appearance than those of 
L. lakedoivnensis, probably because the hairs are 
longer in L.forresti rather than being more 
numerous per scale. In general the tail of both 
species is bicoloured, with the upper part similar in 
colour to the dorsal fur and the lower part paler, 
similar in colour to the ventrum fur. Tails of 
L. lakedoivnensis range from only slightly bicoloured 
to distinctly bicoloured. In general the pes of both 
species are cream to white, with the hair on the feet 
usually the same colour as the ventral hair. 

The pattern of pedal variation in the two species 
is the same although there is some individual 
variation. For example, the lower pad under the 
hallux, (the post hallucal pad) is sometimes missing 
or very indistinct in both species. 

Overall, L.forresti has a more chestnut dorsal 
appearance and stark white ventrally. Dorsally, 
base 50% is dark grey with upper half, buff; 
ventrally, snow white. Tail with almost no bicolour 
and some tails hairier than others. 

The pelage of L. lakedoivnensis varies across its 
geographic range. In the Kimberley, at the Mitchell 
Plateau the dorsum is grey with chestnut tips; the 
ventrum is grey with chestnut tips; the tail is 
bicoloured, pale below, brownish above; and the 
footpads are tiny, sometimes with one indistinct. At 
Mt Percy, the dorsal fur is chestnut tipped (50% to 
33 % of length) over grey; the ventrum is cream; the 
tail is slightly bicoloured, brown above, cream 
below; and the footpads are the same as for the 
Mitchell Plateau population. At Lissadell the dorsal 
fur is chestnut tipped (33 % of length) over grey; 
the ventrum is white to creamy with a little grey; 
the tail is slightly bicoloured to unicoloured; and 
the footpads are the same as for Mitchell Plateau 
specimens. 

In the Pilbara region, the Thevenard Island 
specimens have a dorsum that is either chestnut 
(33%) over brown or all chestnut; the ventrum is 
either brown with the tip cream or is totally cream; 
the tail is unicoloured or weakly bicoloured; the 
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footpads are similar in pattern to other 
L. lakedoivnensis and variable. At Tom Price the 
dorsum is chestnut (33%) over grey but with an 
additional grey tip; the ventrum is grey/cream; the 
tail is bicoloured, dark grey above, paler below. At 
Millstream the dorsum is 25% chestnut over dark 
grey; the ventrum is chestnut/grey or mainly cream 
with some grey; the tail is bicoloured, grey above, 
pale below. At Newman the dorsum is 25% 
chestnut over dark grey; the ventrum is chestnut/ 
grey; the tail is bicoloured as for dorsal and ventral 
fur. At West Angeles the dorsum 25% chestnut over 
grey; the ventrum is chestnut/grey; the tail is 
bicoloured as for the dorsal and ventral fur. 

In the Great Sandy Desert the dorsum is almost 
entirely chestnut; the ventrum almost entirely 
white; the tail is bicoloured, cream below and pale 
brown above; the footpads are very small - the post 
hallucal pad is indistinct. 

In Queensland specimens the dorsum ranges 
between all chestnut to chestnut (50%) over dark 
grey, through to chestnut (50%) over pale grey; the 
ventrum ranges from all white to cream; the tail 
ranges from slightly bicoloured to distinctly 
bicoloured; and the footpads range from a small 
indistinct post hallucal pad to distinct post hallucal 
pad. 

In the Northern Territory the dorsum ranges from 
chestnut (50%) over dark grey with to (50%) 
chestnut over pale grey; the ventrum ranges from 
white (66%) over dark grey to white (66%) over pale 
grey; the tail ranges from slightly bicoloured to 
distinctly bicoloured; and the footpads range from 
small indistinct post hallucal pad to distinct post 
hallucal pad. 

DISCUSSION 

Prior to 1997, only two specimens of Leggadina 
had been collected from the Pilbara region. 
However, in the past five years, various faunal 
surveys by environmental consultants and 
government scientists have seen some thirty 
additional specimens collected and lodged in the 
Western Australian Museum, with other 
individuals marked and released. It appears, 
therefore, that L lakedoivnensis is more widespread 
and common than previously thought, and perhaps 
undergoes periods of both population increase and 
scarcity (Stuart Anstee, pers.com). 

A previous study of genetic variation within 
Leggadina (Baverstock et al. 1976) examined the 
karyotype and variation at sixteen allozyme loci, 
with L. lakedoivnensis represented by only three 
specimens from Queensland. Using a more 
geographically diverse range of specimens and 
many more allozyme loci, this study confirms that 
the level of genetic differentiation between 
L. lakedoivnensis and L.forresti is substantial and falls 
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well beyond the range of intra-specific mammalian 
variation (Nei 1987). Furthermore, the mean 
heterozygosity within each of the populations 
accords well with other mammalian species (Nei 
1987). The single specimen sampled from the 
Thevenard Island population was homozygous at 
all loci; this may reflect genetic drift in an isolate or 
it may be due to sampling error. 

This study indicates that genetic differentiation, 
as measured by 44 nuclear loci, is minimal 
between populations of L lakedownensis across its 
range. Moro et al. (1998) reported mtDNA 
differences between populations of L lakedownensis 
and noted that the Pilbara and Thevenard Island 
populations showed the greatest divergence from 
the other regions. This dichotomy is not supported 
by the allozyme data, which instead group all 
three Western Australian populations together to 
the exclusion of those to the east, and suggest that 
the Queensland population (albeit represented by 
a single specimen) is the most differentiated. The 
Pilbara/Thevenard Island animals are 
considerably larger than all others; indeed, the 
Thevenard Island animals in particular are so large 
that some have suggested that they may represent 
either a distinct species or a geographic outlier of 
L. forresti. Neither of these suggestions are 
supported by any of the three data sets now 
available. Of course, it might still be argued that 
the Pilbara/Thevenard Island populations warrant 
recognition as a distinct subspecies, but the 
overlap in sizes between the Pilbara and Mount 
Percy populations, coupled with the lack of any 
allozyme differentiation in Western Australia, 
leads us to conclude otherwise. 

The skulls of L lakedownensis and L forresti can be 
distinguished using the ratio of M^L/M^L and on 
the degree of elongation of the anterior cingulum 
on in L lakedownensis. Externally, L. forresti can 
easily be distinguished from L lakedownensis by its 
more rufous dorsal and pure white venter. L. forresti 
also has longer ears than any population of 
L. lakedownensis. 

The skulls of L lakedoivnensis morphologically are 
quite variable including features that might 
normally be expected to distinguish murine species 
(e.g. incisor orientation). However, this species has 
the largest geographic distribution of the two, and 
therefore experiences a wide variety of habitats and 
climates. The Thevenard Island population has a 
much larger body size than the other populations. 
Pelage pattern also shows marked differences 
between populations and, together with body size, 
can be used to confidently assign animals into their 
appropriate geographic region. L forresti has a more 
discrete distribution in central Australia, and this 
appears to be reflected in the morphologically more 
homogenous nature of this species. 

According to Reid and Morton (1995), the 
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distribution of L. forresti extends further into 
Queensland and New South Wales than has been 
surveyed herein. Accordingly, we may not have 
fully characterized the genetic and morphological 
diversity present in this species. Based on present 
knowledge of extant populations, the geographic 
distributions of tl'ie two Leggadina species are clearly 
disjunct, permitting species identification in the 
field to be predicted from an animal's geographic 
location. Fortunately, sufficient differences in their 
external morphology are available to make this a 
method of convenience rather than one of need. 

In summary, the morphological and allozyme 
data presented herein demonstrate that the genus 
Leggadina comprises two allopatric species, 
L. lakedownensis and L. forresti. These can be 
diagnosed using both morphological and genetic 
criteria. Of the two, L. lakedownensis displays more 
geographic variation than L. forresti, but this 
variation is nevertheless insufficient to warrant 
subspecific status for any regional populations. 
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Appendix 1 Collection data for specimens used in genetic analyses. Prefixes for registration numbers are SAM - South 
Australian Museum; WAM - Western Australian Museum; NTM - Museum of the Northern Territory; JM 
- Queensland Museum; ABIC - Australian Biological Tissues Collection 


L.forresti 
Northern Territory 

MacDonnell Ranges, 23°26’S, 134°44’E, ABTC24040 
South Australia 

Yelpawaralinwa Water, 27°05'S, 138®4rE, SAM M17585; 
Koonchera Waterhole, 23°43'S, 139°31'E, SAM M17570; 
Dulkaninna - 28°59'S, 138°33’E, SAM M17835, 28°57S, 
138°38’E , SAM M17834; Dalarinna Hill, 29°4rS, 138°41’E, 
SAM M17847; North Olary Plains-32^02'S, 140°54’E, SAM 
M18869, 30°28‘S, 140°57'E, SAM M18876; Todmorden 
Station, 27^22’S, 134°3rE, SAM M19214; Nilpinna 
Station, 28°13'S, 135°25'E, SAM M19251; Pitpalyatjara, 
26°10'S, 129°09’E, SAM M18343-4 


L. lakedozvnensis 
Northern Territory 

Litchfield N.P., 13°29’S, 130°42’E, NTM U4198; Bradshaw 
Stn, 15°23'S, 130‘^09’E, ABTC30310; Gregory Creek, 
15n9'S, 13r08’E NTM U386; Kakadu, 12°40’S, 132°50’E, 
ABTC28413-4 

Western Australia 
Kimberley 

Mitchell Plateau - 14'^52'05‘'S, 125°50'45”E, WAM 
M21786, 14'=*49'15*’S, 125°50'15"E, WAM M21762-3; 
Lissadell-15°23'S, 130°09'E, ABTC30310, 15n9^S, 
130°08'E, NTM U386 

Pilbara 

Millstream - 21°24'49”S, 117°09'43’'E, WAM M47733, 
2n9'48''S, 117°ir20"E, WAM M47720-1, 2r35'S, 
117°04'E, ABTC10610; Mt Tom Price, 22°44’S, 117°47'E, 
WAM M47829; Newman, 22‘^55'S, 118‘^43’E, WAM 
M47912; Thevenard Island, 2r28'S, 115°00'E, ABTC10622 

Queensland 

Paluma, 19°00‘S, 146°0rE, JM10863 
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Appendix 2 Collection Data for specimens used in morphological study. Locality, geographic coordinates, sex (M 
male, F female) and registration number. Prefixes for registration numbers are as in Appendix 1. An 
asterisk beside the Museum number denotes this specimen also used in genetic analysis. 


L.forresti 
Northern Territory 

MacDonald, 22°25‘36"S, 135°07'36"E, f, NTM U2958; 
Oodnadatta, 24°48’00"S, 135‘^22'36’’E ,m, NTM U6352; 
Coorabulka, 23°44’48"S, 140n8'00"E, f, NTM U9688; 
Yuendumu, 22°15'00'’S, 13r48’00"E, 2f, NTM U13628-9, f, 
NTM U13797; Hamilton Downs, 23°3r36"S, 133°16'36"E, 
m, NTM U13632 

Queensland 

Birdsville , 22°23'48"S, 139°38'48"E, 2m, NTM U6350-1 
South Australia 

Innamincka, 27°44'48"S, 140°46’36"E, m, SAM Ml0170; 
Embarka, 27°43'36"S, 140°08'48"E, m, SAM M11476; 
Ucatunna Hill, 27‘’24'02"S, 135°52"49", m, SAM M18514; 
Unducarra Hill, 26°28'32"S, 134°19'53’'E, f, SAM M18545; 
Memory Bore, 26°42’06"S, 135°30'00'*E, f. SAM Ml 8563; 
Four Hills Trig, 28°30'00'’S, 136°30’00"E, m, SAM M18710; 
Abertachna Well, 29‘^43T2”S, 137°20’48''E, m, SAM 
M18808; Emerald Spring, 29'^23'00"S, 137°08’48"E, m, 
SAM M18816; Starvation Dam, 30°27'47"S, 140°56’48"E, 
m, SAM M18876; McIntyre Bore, 30°25’57"S, 140°58‘26"E, 
m, SAM M18882 

Western Australia 

SW Mann Range, 26°18'00"S, 128°52’36''E, 2f, SAM 
M6344-5; Bell Rock, 26®18'00”S, 128°52'36"E, 2f, SAM 
M6347, SAM M6349 

New South Wales 

Fowlers Gap, 31°00’00"S, 141°42'00"E, m, SAM M13633 

L. lakedownettsis 
Western Australia 
Kimberley 

Lissadell - 16°44’00"S, 128°26‘36''E, 4m, WAM M16893, 
WAM16869, WAM M19278-9, f, WAM M16870; 
16°40’00”S, 128°23T3’‘E, 3m, WAM M16872, WAM 
M16874, WAM M16900; 16°44'00"S, 128°26'36"E, m, 
WAM M16893; 16°32'5r'S, 128°33‘30"E, m, WAM 
M18833; 17°37'00"S, 124‘’5219"E, 3m, WAM M16904-6; 
Mt Percy- 17°35’54"S, 124°54'06’'E, 2f, WAM M16919-20; 
17°4r32"S, 124°56’03"E, m, WAM M16944, f, WAM 
M16945; 17°40'23"S, 124®56'03"E, 3m, WAM M16944, 
WAM M16950-1 2f, WAM M16945-6; 17°40'23"S, 
124°56’03’*E, m, WAM M16957, f, WAM M16959; 
17°37'00'’S, 124°52'19"E, 3f, WAM M16952, WAM 
M16962, WAM M16975; Barker Gorge- 17°14’15"S, 
124®38T5'’E, f, WAM M18553, 6m, WAM M18554-7, 
WAM M18559-60; Calder River, 16°19'00"S, 124°58'00''E, 
f, WAM M17494; Mt North , 17°29'46’'S ,124°45'05"E, f. 


WAM M16965; Mitchell Plateau- 14°53‘59’'S, 125°50T0, 
m, WAM M15536, 14°35T5"S, 125°45'30‘'E, 2f, WAM 
M18068-9, 2m, WAM M18070, WAM M18084; 14°47'59"S, 
125°49'00"E, m, WAM M18100; 14"52'05"S, 125°50'45"E, f, 
■^WAM M21526; 14°53'09"S 125‘^48’55"E,m, WAM 
M21623; 14°49T4"S, 125°50'15"E, m, "^WAM M21762, 2f, 
WAM M21752, WAM M21767; 14°50'30"S, 125°50'30’'E ,f, 
*WAM M21786; 14°52T9"S, 125°49T5"E, 3f, 

WAMMM22044, WAM M22090-1; Kununurra- 
15°47'30''S, 128°45'30"E, m, WAM M14857; 15‘’46'00"S, 
128°44'00''E, m, WAM M19644, f, WAM M19645; Bungle 
Bungle, 17‘’15'00"S, 128n8'00"E, f, WAM M32657; 
Mclarty, 19°30'00'’S, 123°30'00"E, m, WAM M19840; Gt 
Sandy Desert- 20‘^24'00’'S, 122“08'00’*E, 3m, WAM M2725, 
WAM M22725-6; 22°33'00"S, 122°23‘00'’E, m, WAM 
M22727; 20°20'59'‘S, 127°26‘00'’E, m, WAM M22728 

Pilbara 

Millstream- 2r35’59"S, 117°04'00''E, f, WAM M27011; 2m, 
■^WAM M47720-21; Tom Price- 22‘^45'00"S, 117°46’00"E, m, 
WAM M47326, 22°49’01"S, 115°46'5TE, f, WAM M47710, 
22°48'48”S, 117°47'20"E, 2f, WAM M47711-2; Chichester, 
2r25’48''S, 117°10’43"E, m, nVAM M47733; Newman, 
22°55'57’S, 118°43’48”E, f, *WAM M47912; Westangeles- 
23ni'5r‘S, 118°45'48‘'E,m, WAM M47672, 23°ir5r’S, 
118'^44'08"E, f, WAM M47673 

Thevenard Island 

Thevenard Island, 2r28'00'’S, 115°00'00"E, 9m, WAM 
M7402, WAM M7404-6, WAM M7408-10, WAM M7412, 
WAM M28921, 9f, WAM M7403, WAM M7411, WAM 
M21296, WAM M34163, WAM M47416-7, WAM 
M47434-5 

Queensland 

Lakeland,15°50'48"S, 144°50’48’’E, 2m, JM1024, JM2874, 
2f, JM1293-4; Mapoon, 11°59'48’'S, 142W00"E, JM4030; 
Batavia, 12°09WS, 141°54'00"E, JM4031; Stonexing, 
12°23’48"S, 142W00"E, f, JM4032; Hillgrove, 19°38'48’'S, 
145°43'12"E, m, JM4954; Hiddenvall, 19°00'00"S, 

146°05’48’’E, f, JM10863; Undurra, 18°13’12"S, 

144°34'12”E, f, JM11573; Rolleston, 24°28T2”S, 

148°38'48"E, f, JM11585; Williams, 14®37'12’‘S, 

144“06‘00"E, f, JM17919 

Northern Territory 

Pine Creek , 13°30'00"S, 13r28T2"E , f, NTM U2445; 
Litchfield, 13W00"S, 130°26'24"E , f, NTM U4198; 
16°04T2"S, 130n5'00"E, m, NTM U4353; Mt Sanford, 
16°35'24"S, 130°20'24’'E, m, NTM U4480; Kakadu- 
13°25'00"S, 132°06’00"E, f, WAM M43288, 13°29’59"S, 
132'^15’00"E, f, WAM M43289, 13‘’30'00''S, 132°35'00"E, f, 
WAM M43290, 13°3r00"S, 132°27'00'’E, m, WAM 
M43291,13°23'59"S, 132n6‘00’'E, m, WAM M43292. 



Appendix 3 Summary statistics for twenty skull and six external variables in eight populations of Leggadina lakedownensis and in Lforresti. (F = female; M = male) 
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